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Abstract 

In this paper, we study the progressive collapse of 3D framed structures made of reinforced concrete 
after the sudden loss of a column. The structures are represented by elasto-plastic Euler Bernoulli beams 
with elongation-rotation failure threshold. We performed simulations using the Discrete Element Method 
considering inelastic collisions between the structural elements. The results show what collapse initiation 
and impact-driven propagation mechanisms are activated in structures with different geometric and 
mechanical features. Namely, we investigate the influence of the cross sectional size and reinforcement a 
and of the plastic capacity /3 of the structural elements. We also study the hnal collapse extent and the 
fragment size distribution and their relation to a, p and to the observed collapse mechanisms. Finally, we 
compare the damage response of structures with symmetric and asymmetric reinforcement in the beams. 
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1 Introduction 


Local damage of buildings can either be due to accidental events like gas explosions, gross design-construction 
errors and malicious terrorist attacks (Levy and Salvadori, 1992), or can be thoroughly planned as part of 


controlled demolition processes with blast. Subsequent cascades of failures can cause large economic and 
human loss when triggered by accidental events and make the difference between effective and dangerously 
ineffective demolitions. 

Research about progressive collapse of buildings proceeded discontinuously since 1970s mostly prompted 
by outstanding and shocking catastrophes. Interest in the subject rose after the Ronan Point partial collapse 
in 1968 due to gas explosion. During the seventies, the fundamental approaches to structural robustness as 
well as many indicators, like the Reserve Strength Ratio (RSR) ( |Maes et al. 2006), were formulated, also 
with regard to off-shore structures that suffered brittle collapses on the North Sea. Renewed attention to 
the problem was given due to terrorist attacks against the Alfred P. Murrah federal building, Oklahoma City 
1995, and the World Trade Center (WTC), New York 2001 ( Val and Val[ 2006). 

Nowadays many codes prescribe alternate paths for the load (Alternate Load Path Method (ALPM) 


(Val and Val 2006)) and high toughness of structural members and their interconnections (Nair 
Nonetheless, these measures are not always sufficient to prevent progressive collapse (Vlassis et al. 


2004). 


2006). 


Moreover, even though the serious damage amplification due to dynamics has been exhaustively pointed out 
in (Pretlove et al. 1991 Marjanishvili and Agnew 2006), static analyses are still used in the context of the 
ALPM. 
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Developing efficient tools to evaluate structural robustness and to prove the effectiveness of measures 
aimed at preventing progressive collapse is therefore an important issue and today several algorithms and 
models are available in literature. A simplified approach to take dynamics and impacts between falling 
elements into account was proposed by (Vlassis 20071 and (Chiaia and Masoero 2008) showed its analogy 
with the variational approach to fracture mechanics. The scheme is based on energy balance, requires only 


static analyses, and was effectively applied to strain hardening (Vlassis et ah, 20061 and softening (Chiaia 


and Masoero 2008) structural elements. 

Analytical ID models were developed after the WTC collapse (]Bazant and Verdure 2007 Cherepanov 


and Esparragoza, 2007 Seffen, 20081. In these models, progressive buckling of the columns is due to the 


impact of the upper floors, considered as an increasing falling mass. Differently, computer simulations permit 
to study more complex 2D and 3D structures. Several key factors that influence robustness of frames have 
already been identified. For instance, we know that the loss of external columns from the facades or the 
corners of a buildings are the most serious scenarios where, according to the ALPM method, one column 
is instantaneously removed (see, e.g. (Kaewkulchai and Williamson 20031). Moreover, it was shown that 
beam-column connections are critical points of failure initiation (Khandelwal et al. 2008) and that catenary 


effects in the floor slabs can remarkably improve robustness (Vlassis et al. 2006). 


Even though the final outcome of progressive collapse depends on the collisions between structural ele¬ 
ments, most of literature focuses on collapse initiation. Collisions are rarely taken into account either detailed 


with Finite Elements (Hartmann et al. 2008 Luccioni et al., 2004), or approximated in the framework of 


Finite Macro-Elements (Kaewkulchai and Williamson 2006 Grierson et al. 2005). Detailed Finite Elements 
are too demanding in terms of computational time for extensive parametric studies on large structures. Dif¬ 
ferently, Finite Macro-Elements are efficient and can be applied to large structures, but they require strong 
approximations to take into account collisions and catenary effects, especially in 3D (e.g. see (Isobe and 


Tsuda, 2004)). 


The lack of experimental results of progressive collapses suggests an approached based on simulations 
whose reliability arises from the basic physics incorporated. The results obtained with such algorithms can 
be used to construct, test and calibrate simpler models. In this work, we use spherical Discrete Elements 
(DE) to simulate the progressive collapse of typical 3D framed structures made of reinforced concrete with 
fixed regular overall geometry (see Sec. [^. The aim is to study the collapse initiation mechanisms due 
to dynamic stress redistribution, and the subsequent damage propagation mechanisms due to collisions 
between the structural elements. Understanding the activated mechanisms, depending on the strength, the 
stiffness and the plastic properties of the structural elements, can help to choose optimal robustness oriented 
design solutions as well as the most appropriate structural reinforcement of existing buildings. We perform 
parametric studies scaling the cross sectional size and reinforcement by the cross sectional scale factor a 
and varying the plastic capacity of the structural elements (see Sec.[^. In this way, we show the expected 
collapse mechanisms and the final consequences of progressive collapse in terms of final collapse extent and 
fragment size distribution for various (a, j3). Finally, in Sec.j^we compare the damage response of structures 
with symmetric and asymmetric reinforcement in the beams. 


2 Simulating Progressive Collapse 


The choice of spherical DE as simulation tool is motivated by several factors (Poschel and Schwager 2005): 


first of all, DE are naturally suitable to deal with dynamic problems since they are based on the direct inte¬ 
gration of Newton’s equations of motion, which makes the algorithm simple and fast. Moreover, geometric 
and material nonlinearities, as well as local ruptures can be easily modelled without remeshing. Momentum 


transmissions due to collisions can be included straightforwardly (see Sec. 2.2). A quite fine mesh is required 
to represent the actual volume of the structure and to reduce the error originating from the fact that instead 
of considering sectional ruptures we instantly remove beam elements that are responsible for the cohesion of 


the system (see Sec. 2.2). Considering sectional ruptures would require remeshing while a more precise rep¬ 


resentation of the volumes could be obtained with polyhedrical DE. For both strategies, the computational 
demand would grow remarkably. Our model has previously been employed to study fragmentation of mate- 

and its applicability to progressive collapse of structures is demonstrated 


rials, e.g. ( 

Carmona et al. 2008 

in (Masoero et al. 

2010 

)• 
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In the simulations, the intact structures are first equilibrated under external service load and gravity. 
If some elements fail during this initial phase, the structure is incapable of carrying the service load and 
the simulation is stopped. Differently, if no elements fail, the local damage induced by an accidental event 
is considered by a sudden removal of a central column of a facade at the first floor (Fig. [^, according to 
the ALPM. The subsequent dynamic stress redistribution can break other elements and trigger widespread 
progressive collapse. The dynamics of the system is followed by means of explicit time integration, using 
a order Gear predictor-corrector scheme that, for the explor ed s et o f pa rameters, is stable with time 
increments lower than 10“® seconds. In the following subsections 
of the studied structures and employed model. 


2.1 


and 2.2 we give a detailed description 
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Figure 1: The studied framed structure with a) small and b) large a. c) Cross sections of the structural 
elements and arrangement of the reinforcement. 


2.1 Construction of the framed structures 

For comparative reasons we limit ourselves to study typical regular 3D frames formed by 4x4x4 identical 
square cuboid cells with L=4m and iJ=3m (see Fig. [^a). The structures are made of columns along the 
vertical Z direction, clamped to the ground a,t Z = Om and connected at each storey by principal beams in 
X and Y direction. Thin slabs spanning between the principal beams form the floors while the presence of 
walls is not considered. The geometry of the cross sections of the structural elements is displayed in Fig.j^c, 
where the subscripts c, b and s denote columns, beams, and slabs. We set the height of the cross sections 
proportional to the length of the structural element with the coefficients A, namely /ic=H/10, /it,=L/I0 
and /is=L/50, and we scale each /i by a dimensionless cross sectional scale factor a. a is identical for all 
elements and enlarges their cross sections making the structure stiffer and stronger. The base edges bb, 
be of the cross section are proportional to the heights hb, be with the aspect ratio coefficients Sb=2/3 and 
(5c=l- Consequently, the area, the sectional inertia with respect to the f principal direction and the torsional 
inertia of the cross sections of beams, columns and portions of floor slabs can be easily computed. We 
represent a structure made of reinforced concrete (RC) with Young’s modulus Ee and shear modulus Ge (see 
Appendix, Tablej^. The reinforcement is symmetrically distributed (see Fig.j^c) and its area is proportional 
to that of the cross section by the factors ps_c=l-78%, 6=0. 58%, and ps^s=1.26%. The structure carries its 
own weight G, the service external dead load Gd=285kg/m^ given by non structural elements like pavement, 
plaster and internal walls, and the service live load (5=667N/m^. Gd and Q are considered uniformly applied 
to the upper faces of the floors. 
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2.2 Model Description 

We represent the columns and the beams by meshes with ric and Ub Euler-Bernoulli (EB) beam elements 
respectively. The floor slabs are considered by a grid of Ub x Ub EB elements of length Lg (see Fig. |b) that 
define slab portions of size Lg x Ls x hg (see Fig. [^c). To represent the volume of the structure, we set Ub 
and ric to obtain « 1, where /i® and L® are the height of the cross section along r] and the length of the 
generic EB element. The cross sections of the EB beam elements are set according to Sec. |2.1| The error on 


Irj introduced by the simplifying hypothesis Irj = i 
the slabs in the horizontal plane is not relevant. 


a) 



acceptable because the bending of the beams and of 

b). 



Figure 2: a) Spheres and b) EB elements in the modelled structure with a = 1.43. The light grey area marks 
the initially removed column. Generic Euler-Bernoulli element in c) undeformed and d) deformed state. 


We represent the structural volume by spheres surrounding each node (see Fig. [^a). The diameter of 
the sphere is equal to 90% of the length of the shortest EB element connected to node k. The mass Mk 
is obtained summing the contributions from the mass Meb oi the m EB elements connected to node k and 
from an extra mass Mex,k given by the external dead load Gd' 

m ^ 

Mk = ■ ( 1 ) 

Note that in the DE algorithm the EB beam elements do not have a mass since Meb is concentrated in 
the spheres. The rotational inertia of a sphere is computed considering Mk to be uniformly distributed. 

2.2.1 Euler-Bernoulli beam elements 

The EB elements determine the interactions between pairs of nodes, associating their relative rotations and 
displacements to forces and moments acting on them. In the following, we describe the linear elastic - 
perfectly plastic constitutive behavior and the failure rules of the EB elements. 

Elastic regime In the linear elastic regime, we use the force-displacement relations described in detail in 
2008[ ), thus taking into account the geometric nonlinearities due to large displacements and 
neglecting shear deformability. The rotations yi are defined starting from the deformed state of the generic 
EB element (see Fig. j^d). Namely, the bending rotations (prj around the rj principal axes align the ^o, 
Cl axes with the line connecting nodes 0 and 1, while rotating the Cii Vi axes around C by makes them 


(Carmona et al. 
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parallel to ^q, tiq: 
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(2) 


At a given time step, we compute the axial strain e = AL’^/Lq and the rotations (p for every EB element 
starting from the position of the spheres and at the orientation of the C, i?, C axes frozen to nodes 0 and 1. 
The forces and moments at nodes 0 and 1 come from the EB beam theory: 


Bo,3 
Bi,3 
Tz,3 


N 


EJn EJf ff 

(4(po.i + 2ipi,j) = 

re 

j-'n 


E, 


+ 2(/3oj) = 
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^0 


Eo,j + B 


ij 


±EcA^e , 

,GJ! 


(3) 

(4) 

(5) 

( 6 ) 
(7) 


Ec is the Young’s modulus of concrete (see Table [^. Bi j and (p^Y denote the bending moment and the 
effective rotation around the j = CtV axis at node z = 0,1, Tij is the shear force along the j axis at node 
TV is the normal force and Mi^ is the torsion. Damping inside the beams is considered adding forces and 
moments at nodes 0 and 1 that are proportional to the elastic part of the velocities e, po,^, Po.rj, 
ip(^ by the coefficients in Table but with opposing direction 


Poschel and Schwager |2005 ). 


Plastic regime Progressive collapse of structures involves many irreversible processes in the elements and 
plastic energy dissipation can determine robust or vulnerable responses to damage. To consider plasticity in 
the EB elements, we make the simplifying assumption that axial and bending plasticization are uncoupled. 
This choice is justified by our intention to keep the model as simple as possible, leaving refinements to 
further works. Furthermore, we don’t consider plasticization due to shear or torsion in the RC because 
plastic dissipations associated to them are generally small. Under these hypotheses, the EB elements enter 
the perfectly plastic regime in axial direction or in bending at node i and around ^ or yy independently if one 
or more of these conditions is satisfied: 


N>Ny or N < -|TV,^| , (8) 

\B.,3\>By . (9) 

We set the tension and compression yield thresholds = p^A^^fy, N]! = A^^fc neglecting the contribution 
of concrete in tension and steel in compression. The bending yield threshold is evaluated referring to the 
^ axis and neglecting the contribution of concrete. If is the fraction of reinforcement in tension, i.e. | for 
the columns and ^ for the beams and the slabs (see Fig.[^c), we obtain: 

By = tipiA^fyh^ . ( 10 ) 


We also add a further contribution to By to consider that compressions N < 0 increase By by 

reducing the area of concrete in tension during bending. We set AM^ assuming that it is carried by the 
reinforcement alone and that it compensates the strain Sg in the reinforcement in tension produced by N, 
i.e. es{N) = es{AMy): 


N AMy 

A^Ec ~ EgPlA^h^Eg 


( 11 ) 


Within the employed direct time integration scheme, we check yielding in terms of strain and rotations 
instead of forces and moments. Thus, we adopt elongation e^, shortening and bending p^ff,y yield 
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thresholds that satisfy Eqs. E-E in the equality form, when inserted into Eqs. [^-[^ The expressions for the 
and are summarized in Table where also the additional term to from Eq. 11 is shown. 

In the perfectly plastic regime, we consider the axial strain e and rotations (po, (pi to result from the sum 
of an elastic and a plastic contribution. Fig. shows how axial plasticization is implemented at the generic 
time step ti. History dependence is considered accumulating the plastic strain in time. 
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AeP' = 0 


E®' = E^ 

AeP' = £ - e'' 




£P'(t) = EP'(t )+A£P' 


From DE step 


e(tj = £(Fi)-£P'(t.) 


Figure 3: Operations to check axial plasticization and to compute the axial plastic flow. 


The linear distribution of bending along the EB element makes the description of the rotational plastic 
regime more complicated. First, we compare ^p^J^■, ^p^x^ obtained from the integration of Newton’s equation 
with ip^ff-y to define whether plasticization occurs only at one node or at both. If for instance only node 
0 enters the plastic regime, then we set (po = ipQ + , where (pQ , put into Eq. satisfies Eq. in the 

equality form. Differently, if both node 0 and 1 plasticize, either pQ and pi must be separated into elastic 
and plastic part, so that Eq. |^and Eq. [^in the equality form return a linear system of two equations in the 
two unknowns Pq^ and pf. As for axial plasticization, the plastic parts Pq^ and Pi^ must be subtracted from 
the total Po and pi at the next time step, and the plastic rotations must be cumulated in time. 


Element failure If the strain in a cross section of an EB elements is too high, the element fails and is 
instantly removed from the system. In the following, we will scale the plastic capacity of the EB elements 
by a parameter /3 >0. If /I = 0 the elements break when a combination of the elongation and the effective 
rotations at nodes z = 0,1 is large with respect to the yielding thresholds evaluated in uncoupled conditions: 
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( 12 ) 

(13) 


(14) 

(15) 


with ultimate threshold values of elongation e*^, shortening e*^ and rotation p^^ estimated in uncoupled 
conditions. Failure due to shear and torsion is neglected because the shear reinforcement is thought to be 
designed according to the capacity design approach, that avoids the occurrence of these brittle mechanisms 
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before bending or axial strain failure. We assume and equal to the ultimate tensile strain of steel 
and compressive strain of concrete respectively (see Table [^. is estimated considering a state of uniform 

bending, and thus uniform curvature x, in the EB element. Under these hypotheses, the rotation between 
the edges of an EB element is x-b® and gives a strain in the steel bars at ry = /i®/2 equal to = x ■ ^^/2. We 
thus obtain setting Eg equal to the ultimate strain of the steel, i.e. Eg = Eu,s- The expressions in Table 
[^consider also that the adopted meshing rule assures« 1. 


2.2.2 Hertzian contact between the spheres 

Progressive failure of EB elements can lead to the free fall of structural elements. We model inter-spheres 
collisions by a Hertzian potential ( ]P6schel and Schwager 2005 Carmona et al. 2008) that generates a 
conservative repulsive force between partially overlapping spheres. This force is directed along the line 
connecting the centres of mass of the colliding spheres and is proportional to the overlapping volume by a 
stiffness parameter V. A similar force is generated when a sphere crosses the Z — 0 plane that represents the 
ground. Impacts dissipates energy due to local fragmentation and to sliding and rolling friction. Thus we 
introduce forces that are proportional and opposed to the normal, the tangential and the rotational relative 
velocities of the overlapping spheres with the damping coefficients summarized in Table In tangential 
direction, either Coulomb’s or dynamic sliding friction are considered. 


3 Parametric studies of Progressive Collapse 


The plastic capacity of the structural elements is a key factor of progressive collapse. Plasticity determines 
the subsequent mechanisms of damage propagation based on collisions and the final extent of collapse. We 
employ the presented model to perform parametric studies on the structural cross sectional scale factor a 


(see Sec. 2.1) and on the plastic capacity (3 (see Eqs. 14 and 15). We choose a as parameter because a 
structure with given /3 can be robust or extremely vulnerable to a column removal, depending on the cross 
sectional size of the elements. In this way, we can see the effect of plasticity in structures that exhibit different 
responses to the initial damage, ranging from no collapse to catastrophic collapse. In Sec. |3.1| we describe 
the observed global and local primary mechanisms that can trigger progressive collapse and the subsequent 
collisions-driven mechanisms. In Sec. |3.2| we show the results of the parametric studies. We especially point 
out what collapse mechanisms occur depending on a and /? and what are the consequences in terms of final 
extent of the collapse. Finally, in Sec. |3.3| we show how a and [3 effect the fragment size distribution of the 
rubble. 


3.1 Collapse Mechanisms 

The initial column removal can trigger three different primary collapse mechanisms that start progressive 
collapse, two of which are global and lead to total collapse: 

1. the first one is caused by elastic waves inside the floors that can separate even distant slabs from the 
beams (see Fig. |^a). 

2. The other global mechanism separates less slabs from the beams but the floors are progressively punched 
by the columns (see Fig. ib). Most probably, this mechanism would turn into progressive buckling of 
the columns if they were less reinforced. 

3. The local primary collapse mechanism is characterized by a crack propagating from point A (Fig. |4}c) 
and disconnecting the neighboring floor slabs from the beams (Fig. |^d). 

To explain why rupture occurs at point A instead of point B or C we consider the schematic repre¬ 
sentation in Fig. I^e. After the column removal, the cross sections Ba and Bq, being at two sides of 
the same node B, show the same vertical displacement. In Fig. |^e, the load on qi on beam A — Ba 
is greater than the load (72 on beam C — Be because the area that can transfer load to the former is 
larger. Therefore, considering that the torsional stiffness of C — Be can be represented by a torsional 
spring that reduces the bending moment in Ba and increases that in A, the maximum static bending 
moment is located in A. 
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Figure 4: Global collapse due to a) catastrophic wave (a = 0.51, /3 = 0.2) and b) progressive punching 
(a = 0.54, P = 0.8). Local starting collapse (a = 0.56, /3 = 0.2): a) detail of the first failure area, where the 
light gray area marks the initially removed column and the arrows show the direction of crack propagation; 
b) first stages of the local progressive collapse, e) Approximate static schemes of the elements where the 
starting damage propagation occurs; the springs in B represent the stiffness of the perpendicular beams. 


When a local primary collapse mechanism is triggered, the portion of structure above the removed column 
undergoes free fall and collides with the floor slabs below. This hammer effect (Fig. generates elastic 
waves that can damage the neighbouring floor slabs. Rarely does it have catastrophic consequences by itself. 




Figure 5: a) Hammer and b) subsequent elastic wave at the first storey {a = 0.54, /3 = 0.2). Free fall of the 
floor slabs, c) almost horizontal for brittle structures {a = 0.59, /3 = 0.2) and d) tilted for plastic structures 
{a = 0.57, P = 0.8). In brittle structures e) the slabs stacking on the ground push less against the column 
than in f) plastic structures 


Other sources of damage transmission are the lateral impacts between falling rubble and still intact 
portions of the structure. This drag effect can destroy the perimeter beams of neighbouring floor slabs, 
which can eventually collapse, but usually is not able to cause a widespread propagation of damage by itself. 
Finally, a severe secondary mechanism is due to the forces exerted by the rubble stacking on the ground, 
which can cut the neighbouring columns at the base {base — cutting, see Fig. ic-f). 
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3.2 Phase diagram of final states 

The collapse mechanisms described in Sec. |3.1| can occur depending on the cross sectional size and reinforce¬ 
ment of the structural elements, determined by the cross sectional scale factor a, and on the available plastic 
capacity j3. Fig. [^a shows the response of the studied frames to the applied local damage. 




Pi/' [rad] 


0.08 


- No total collapse above 

—— No partial coll, below 
• ■ - - No base cutting above 
Sim result: total coll, 
o Sim result: partial coll. 



0.0 0.2 0.4 0.6 

Plastic capacity p 


Figure 6: Response to the initial damage in terms of final collapse. Above the bold dotted line there is not 
base cutting. Above the bold dashed line the eventual primary collapse mechanism is local. Th e mark ers 
denote single realizations. is the maximum plastic rotation in uncoupled condition (cf. Eqs. 


14 


15) 


For (a, P) pairs below the auj{P) curve, the intact structures experience static collapse before the 
damage. Such weak frames are not supposed to exist since they can not carry the service load. Structures 
with (a, P) within the dashed area just above the otuj{P) curve completely collapse after the column removal 
triggers a global primary mechanism. This (a, P) region is narrow and vanishes when P < 0.2 since brittle 
failures induce compartmentalisation, abruptly interrupting the dynamic stress flow. In this way, either the 
propagation of waves (dominating when 0.2 < P < 0.6) and the global failure of the storeys due to progressive 
punching (/3 > 0.6) are avoided. Structures with (a, P) above the dashed region and below the Q!c,n(/3) exhibit 
local primary collapse, but can still experience total collapse prompted by the base — cutting mechanism (see 
Fig. I^b). If P is small, the probability that a sequence of base — cutting provokes total collapse is low, due 
to the almost null tilt of the floor slabs during the free fall, the high degree of fragmentation after collisions 
and the larger a values. On the contrary, when P is large, sequences of base — cutting are frequent because 
the slabs tilt while falling and stack on the ground without considerably fragmenting. Consequently, massive 
slabs lean against the thin (small a) base columns at a remarkable height, and cut them (Fig. [^c-e). 

Partial collapse can occur when /? < 0.8 and a is sufficiently large (see Fig. §b). In this case all the 
collision mechanisms of Sec. |3.1| determine the final extent of the collapse. Overlapping between the partial 
and total collapse regions in Fig. is due to the fact that the consequences of collision driven mechanism 
are considerably variable. Structures with {a, P) pairs above the ac^oiP) curve are perfectly robust, i.e. 
they do not suffer any further failure after the column removal. Finally, when a > ayj the intact structure 
does not plasticize in the pre-dmage stage, i.e. during the static application of the service load. Since this is 
a common requirement for buildings in service conditions, structures with realistic size of the elements are 
located in the a > ayj region. 

The beneficial effect of plasticity is evident from the fact that ac^oiP) and auj{P) decrease with p (see 
Fig. 1^ a). These curves sensibly decrease in the range 0 < /? < 0.2, which means that the provision of a 
minimal plastic capacity is sufficient to remarkably improve the structural response, assuring complete safety 
to structures with a > otyj. This result is a consequence of the symmetric distribution of reinforcement 
inside the beams and the floor slabs, that makes their bending behavior qualitatively similar to that of steel 
elements. Steel structures are actually likely to sustain one column removal without damage propagation. 
In Sec. 1^ we show that RC structures with realistic cross sectional size of the elements and asymmetric 
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Figure 7: a) Final entity of the collapse for {a, j5) pairs between auj{P) and ac,D{l3)- Note that for each 
/3, values of a larger or smaller than those of the plotted points give respectively DR = 0 and DR = 1. b 
Fragment mass distribution for different plastic capacities /3. p{m) is the fraction of fragments with mass m 
normalised to the total mass of the structure Mtot- 


reinforcement distribution would experience partial collapse. Differently, when /3 > 0.4, ac,DiP) and Oc 
do not decrease much anymore (see Fig.j^a). This means that the initiation of progressive collapse is a local 
phenomenon associated with relatively small plastic stress redistributions, as already observed for 2D steel 
frames by (|Khandelwal et al. 2008). 

Fig. [7la shows a quantitative measure of the final collapse extent for structures with different plastic 
capacity /?. In particular, we compute the demolition ratio DR, i.e. the fraction of lost living space at the 
end of the collapse. From the figure it can be immediately seen that structures with high plastic capacity 
/3 undergo progressive collapse only if they have thin elements, i.e. small a, but if the collapse is triggered 
it will affect the whole system. The large variability of DR is due to the collision-driven mechanisms and, 
in particular, the largest values of DR associated with partial collapse, i.e. when DR ^ 0 and DR ^ 1, 
indicate the occurrence of base — cuttings. 


3.3 Fragment Size Distribution 

The size of the fragments produced after a structural collapse is interesting for controlled demolitions, where 
large fragments require further effort to relocate them. The probability density distribution p of the mass of 
the fragments m normalised by the total mass of the structure Mtot is shown in Fig. 0b. We observed that a 
does not influence p{m). Differently, as the plastic capacity of the elements /3 grows the power law regression 
lines in figure shift towards larger sizes of the fragments while their exponent, between -1.2 and -1.4, does 
not seem to change according to a trend. Note that this exponent is in the range of the one reported for 


shell fragmentation of 1.35 (Wittel et al. 2004). When (3 is small p{m) has large dispersion and most of the 
fragments are represented by single spheres completely disconnected from the others. Calling Mi the sum of 
the masses of the fragments made of one sphere and Mp the total mass of the fragments, Mi/Mp is larger 
than 0.7 when /3 = 0. This denotes a finite size effect, i.e. the single sphere is larger than the characteristic 
size of the fragments. On the contrary, when /3 is high p{m) is less dispersed and Mi/Mp < 0.25 implies 
that the fragment size distribution is better caught. 


4 The effect of asymmetric reinforcement 


The results in Sec. 3.2 refer to RC structures with symmetrically reinforced elements. Nevertheless, the rein¬ 
forcement inside the beams and the floor slabs of real structures is mostly concentrated in the regions under 
tension in service conditions. Fig. 0a shows the symmetric and the more realistic asymmetric reinforcement 
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Figure 8: 3D representation of the beams involved in the local primary collapse mechanism in Fig. @c,d) 
Symmetric and asymmetric reinforcement arrangement. Static deformed condition of the beams, bending 
moment in ^ and regions in tension b) before and c) after the column removal. 


arrangement inside the beams that are involved in the local primary collapse mechanism of Fig. |^c,d. The 
column loss in Fig. |^c produces the inversion of the bending moment in Ba and thus, in case of asymmetric 
reinforcement, the reinforcement under tension in Ba passes from Ag before the damage (see Fig. j^b) 
to As after the damage (see Fig. |^c). Therefore, if is small, section Ba is likely to fail before section 
A, differently from the observations in case of symmetric reinforcement (see Sec. 3.1). In the following, de¬ 
pending on 4', we estimate the ac,_D(/3) curves that separates the collapse and robustness regions in Fig.j^a 
for frames with asymmetric reinforcement. 



Figure 9: a) Load scheme ofthe the A — B beam (cf. Eq. 17). b) Transition to no collapse after the column 
removal for different ratios of asymmetry ip of the reinforcement in the beams. 


Failure in Ba occurs if: 

Bb, = ^-^= DYN{P) ■ Bl^ , (16) 

where Bbj, is the static bending moment in Ba after the damage (see Fig.|^c) and DYN{I3) is an amplifying 
factor that considers dynamics and plastic capacity. DYN ranges from ^ for linear elastic - perfectly brittle 
structures (/3 = 0) to 1 for perfectly plastic structures with /3 = oo. g is the load per unit length on the 
beam (see Fig. |^a): 

9 = [Q + (Gd + AsaL7c)g]ro(/3)L , (17) 

Bg^ = is the yield bending moment in Ba, with By oc because it is proportional to the area of the 
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Table 1: DYN/w factors of Eqs. 18 and 19 


p 

0.0 

0.2 

0.4 

0.6 

0.8 

DYN/w 

0.85 

5.00 

6.60 

8.20 

9.80 


cross section times its height (see Eq. 101. Thus, Eq. 16 can be rewritten as. 


[Q + {Gd + \saL^c)g\w{P)^ = DYN ■ msPs,b>?b^bfyO? 


( 18 ) 


For a given 'h, solving Eq. [^in a for different /? permits to trace the Ocn {P) curves. Nevertheless, in Eq. 18 
the parameter DYN/w{P) is still unknown. 

In order to assign values to DYN/vj{P), we repeat the previous argument for structures with symmetric 
reinforcement, for which we now the acj^{P) values at some P obtained from our simulations (see Fig. j^a). 
In case of symmetric reinforcement, the A — Ba beam in Fig. [^c fails in A, where the bending moment in 
qLp/i. Since in this case 'k = 1, Eq. 18 turns into: 


[Q + {Gd + AsaL7c)5]ro(/3)^ = DYN ■ tsPs^^XlSyfya^ 


(19) 


Solving Eq. 19 starting from the ac,^{P) values in Fig. |^a, we obtain the DYN/w{P) values in Table 

Inserting the DYN/w{P) values from Tablej^into Eq.[^ we obtain the curves in Fig.j^b. These curves 
show that collapse can initiate also in well designed RC structures, i.e. with a > ayj. Moreover, we argue 
that the final extent of the collapse should be partial. Note that the q:„,/(/ 3) curve in Fig. does not change 
with ip since it only depends on the reinforcement in tension in the intact structure. Also, the bold dotted 
line in Fig. [^b regarding the base-cutting phenomenon does not depend on ip but only on the reinforcement 
inside the columns, that is usually symmetric. Therefore, when ip < 1 an expansion of the partial collapse 
region towards higher values of a is expected. 


5 Conclusions and Outlook 

Progressive collapse of framed structures after local damage consist in an initial triggering and a subsequent 
damage propagation. If the initial damage is small, like the studied column removal, collapse initiation 
is generally a local phenomenon affecting the surroundings of the initially damaged area. Global primary 
mechanisms can occur only in thin structures with enough plastic capacity to avoid the compartmentalisation 
effect produced by brittle ruptures, i.e. P > 0.2. Nevertheless, if the starting damage is more serious 
than a single column removal, global primary mechanisms can be expected also for larger and more brittle 
structures. In case of multiple column removal, progressive crushing or buckling of the columns is a possible 
global primary collapse mechanism that was not observed in this context. We showed that structures with 
minimal plastic capacity and symmetrically reinforced beams are robust towards a single column removal. 
On the contrary, frames with asymmetric reinforcement would experience partial collapse even if made of 
elements with large plastic capacity. 

If a local primary collapse mechanism is triggered, the final extent of the collapse depends on secondary 
mechanisms driven by collisions between the structural elements. We showed that damage can not widely 
propagate in structures with small plastic capacity since brittle failures compartmentalise the system. Dif¬ 
ferently, structures with large plastic capacity tend to collapse entirely after a sequence of base cutting. 
This result is a consequence of the fact that, in presence of large plastic capacity, we studied frames with 
thin structural elements. In facts, columns with larger and thus more realistic cross sectional size and 
reinforcement would not fail because of base cutting. 

Collision-driven mechanisms also determine the outcome of the fragmentation process. We showed that 
the fragment mass distribution does not depend on the strength and stiffness of the structural elements. It 
is namely influenced by the plastic capacity of the elements. In structures with large plastic capacity the 
fragments are more massive and represent an extra cost in controlled demolitions processes. 
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In the present paper, for clearer interpretability of the results, we limited ourselves to simple geometry, 
collisions, and constitutive models. Implementing more sophisticated collision models, e.g. using polyhedral 
discrete elements, or enabling discrete elements to fragment (see e.g. (Pdschel and Schwager 20051), as well 
as rate effects are future challenges. We also neglected shear failures, since the structural elements were 
sufficiently small and slender, but this hypothesis should be removed to deal with structures made of large 
elements. These models should also be refined if the aim is to simulate in detail the collapse of specific real 
buildings. Nevertheless, already at the present state, several interesting studies can be conducted to analyze 
the response to earthquakes and to investigate the influence of material disorder, geometric uncertainties, 
overall geometry, and structural connections. Experimental validation remains problematic in the field of 
progressive collapse, because of difficulty in monitoring collapse of complex buildings, and due to problems 
concerning repeatability of experiments. Further discussion about monitoring collapse of buildings for model 
validation, especially concerning demolitions, can be found in (]Bazant and Verdure 20071. 
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Tables of parameters 


Table 2: Loads and mechanical properties of the materials and of the frames. 


Parameter 

Symbol 

U nits 

Value 

Properties of the concrete 

Specific weight 

IRC 

kg/m^ 

2500 

Young modulus 

Ec 

N/m^ 

30-10® 

Shear modulus 

Gc 

N/m^ 

15 • 10® 

Compressive yield stress 

fc 

N/m^ 

20-10® 

Ultimate shortening 


- 

0.0035 

Properties of the steel 

Young’s modulus 

Es 

N/m^ 

200-10® 

Yield stress 

fy 

N/m^ 

440-10® 

Yield strain 

^y,s 

- 

0.0022 

Ultimate strain 


- 

0.05 


Table 3: Damping, yielding and failure parameters of the Euler-Bernoulli elements. 


Parameter 

Symbol 

Units 

Value 

Damping coef ficients 

Elongation 

IL 

Ns/m 

100 

Torsion 

It 

Nms 

1 

Bending 

IB 

Nms 

10 

Axial yielding 

Elongation 

ey 

- 

Ny/iE,A^) 

Compression 

£y 

- 

N^liE^A^) 

Effective rotation yielding 

Columns 


rad 

ByLy{EJI)+ 




+£9L^PsE,/{2h^E,) 

Beams & slabs 

l^eff,y 

rad 

ByLy{EJi)+ 




+£6L^PsEs/{h^E,) 

Ultimate thresholds 

Elongation 


- 


Compression 

~th 

- 

^U,C 

Rotation 

^th 

rad 





« 2£u,s =0.1 
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Table 4: Contact parameters. 


Parameter 

Symbol 

Units 

Value 

Overlap stiffness 

Sphere-sphere 

Y 

N/m^ 

w 

Sphere-ground 

Y9 

N/m3 

5-10^ 

Normal damping coef ficients 

Sphere-sphere 

7n 

Ns/m 

5-10^ 

Sphere-ground 

'■yQ 

In 

Ns/m 

10^ 

Sphere — sphere tangential damping coefficients 

Coulomb friction 

h 

Ns/m 

5-10^ 

Dynamic friction 

It 

Ns/m 

5 • 10^ 

Rolling friction 

Iw 

Nms 

5- IQi 

Sphere — ground tangential damping coef ficients 

Coulomb friction 

p9 

Ns/m 

10^ 

Dynamic friction 

i! 

Ns/m 

104 

Rolling friction 

'■yd 
/ W 

Nms 

102 
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